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Abstract 

It is often the case that scalar fields are produced in the early Universe in the 
form of coherent oscillation. These scalar fields may have huge abundances and 
affect the evolution of the Universe. In particular, if the lifetime is long enough, 
they may cause cosmological disasters. We revisit the issue of coherent oscillation 
of the scalar field when it couples with another oscillating scalar field, and find a 
situation that the abundance, or the amplitude of the oscillation, is significantly 
reduced by a variant type of the adiabatic suppression mechanism. As a concrete 
example, it is applied to the saxion, a flat direction in the supersymmetric axion 
model, and we show that the cosmological saxion problem is solved in a particular 
setup. 
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1 Introduction 

In the theory beyond the standard model such as supersymmetry (SUSY), supergravity 
or string theory, various scalar fields appear, which are often called moduli, and they may 
give significant effects on the cosmological scenario [1]. Moduli are in general stabilized 
at the high energy minimum in the early Universe, which are displaced from the low en- 
ergy minimum. At some epoch, the moduli start the coherent oscillation around the true 
minimum with large initial amplitude and hence the Universe may be dominated by the 
moduli. Some of them may acquire the mass only from the SUSY breaking effect and it 
is of the order of the gravitino mass [2] . In addition, they may have interactions with the 
standard model particles typically suppressed by the Planck scale or some high energy 
scales. Hence their lifetimes are long enough and their decays may significantly affect the 
cosmology. The abundance of late-decaying moduli is strongly constrained by observa- 
tions such as the light element abundances synthesized through big bang nucleosynthesis 
(BBN), the distortion of the cosmic microwave background (CMB) spectrum, or diffuse 
X(7) ray background [3, 4]. If the modulus lifetime is longer than the present age of the 
Universe, it contributes to the cold dark matter (CDM) density, and its abundance is 
constrained so as not to exceed the present matter density. 

Aside from the moduli, there may be many singlet scalar fields which more or less have 
similar properties to the moduli. We focus on such a scalar field (f) whose lifetime is so 
long that its decay might affect the cosmology after BBN. Generally, the energy density 
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to entropy ratio of is calculated as 



^ = -T R ( AV ~ 1 x 10 4 Gevf Tfl 1 f AV, (i) 
s 8 R \Mp) \10 5 GeV J \M P J ' 1 J 

for > Tj, where is the mass of 0, T/ is the inflaton decay rate Tr is the reheating 
temperature and fa is the initial amplitude of 0. For < Fj, we have 

^ = I T / r AV~ixlQ»Gev/ r -^-V 1/ V m0 V'Y-^-V ( 2 ) 
s 8 osc VMp; VlOOy 1 VlOOGevJ ^MpJ ' l J 

where T osc represents the temperature when starts to oscillate. Comparing the above 
estimates to the observational constraints, causes cosmological difficulties unless its 
abundance is somehow suppressed. A late-time entropy production process may act as an 
efficient dilution mechanism for such coherent oscillations, but it simultaneously dilutes 
the preexisting baryon asymmetry. Although some elaborate baryogenesis mechanisms 
may work even under the entropy production, we here consider another possibility. 

In this respect, Linde [5] proposed an interesting scenario that the modulus abundance 
is significantly suppressed without any entropy production process. The idea is that if 
the modulus obtains a mass squared of ~ c 2 H 2 where H is the Hubble parameter and c > 
O(10), the modulus adiabatically follows the temporal minimum of the potential toward 
the true minimum. As a result, the amplitude of the coherent oscillation is exponentially 
suppressed. This mechanism has been recently investigated in detail in Ref. [6] where it 
was shown that the modulus oscillation is generically induced at the end of inflation to a 
small but non-negligible amount. Still, the adiabatic suppression mechanism works as an 
efficient way to significantly reduce the modulus abundance for solving the cosmological 
moduli problem [7]. 

The essential point for the adiabatic suppression is that the scalar obtains a effective 
mass larger than the Hubble parameter, which gradually disappears after the Hubble 
friction becomes inefficient. Thus thermal mass for the scalar may also cause an adiabatic 
suppression [8]. In this paper, we propose a variant method to suppress the amplitude 
of the coherent oscillation of the scalar field, although it may not be the moduli in an 
original sense. 

Let us consider a scalar field 0, whose abundance we want to know, and suppose that 
its true minimum is displaced from the initial position determined during/ after inflation. 
In this setup, it is not hard to imagine that the coherent oscillation of is induced 
when the Hubble parameter becomes equal to the mass of and the resulting abundance 
is estimated by (1) and (2). The situation significantly changes if the scalar field 
couples to another scalar field S, which has a large field value and is oscillating. The 
coupling induces an effective mass of 0, and it can be large enough to overcome the 
Hubble parameter. Hence the scalar field follows the time-dependent minimum of the 
potential adiabatically. As a result, the amplitude of the coherent oscillation of is 
continuously damped by the expansion of the Universe and the energy density of is 
highly suppressed compared with the one without the coupling. This phenomenon was 
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already noted in Ref. [9] in a context of SUSY axion model. Here we wish to generalize 
and reformulate the arguments there. 

This paper is organized as follows. In Sec. 2, our proposal to suppress the scalar 
field abundance is discussed by using a simple toy model. Essential ingredients are all 
contained in this simple toy model. An application of our mechanism to the SUSY axion 
model is discussed in Sec. 3. Sec. 4 is devoted to conclusions. 



2 Adiabatic suppression mechanism in a simple toy 
model 

In this section, we illustrate the mechanism to suppress the oscillation amplitude of the 
scalar field. The essence of the idea is that if the mass of the modulus field is much larger 
than the Hubble parameter, the field oscillation follows its time-dependent minimum of 
the potential adiabatically. This is because the potential around the temporal minimum 
is steep enough to overcome the Hubble friction. Then the resultant amplitude of the 
oscillation becomes very small compared with the general case considered in the previous 
section. Here we introduce an another scalar field, which supplies the large mass to the 
scalar field. To see the suppression of the amplitude, we follow the scalar field dynamics 
within a simple toy model. 



2.1 Toy model 

Let us consider a simple toy model, where two gauge singlet real scalar fields are intro- 
duced. One is the field <fi whose abundance is what we focus on, and the other is the 
heavy field S which supplies a large mass to <p. The scalar potential is given by 

^=^(0-0 o ) 2 + ^ 2 + iAW, (3) 

where and mj are the masses of <j) and S, respectively and A is a coupling constant 
which is assumed to be positive for simplicity. This is the most general form respecting 
the Z2 symmetry under which S transforms as S — > — S, up to the quadratic terms in 
<p and S. The true minimum of the potential lies at S — and <fi = <p$. In the early 
Universe, however, both <p and S may be displaced from the minimum. In particular, if 
S has a large amplitude, the temporal minimum of 4> is calculated as 

2 

771 

<*> = sj+k*- (4) 

This is time-dependent and estimated as ((f)) ~ for \XS\ ^> and ((f)) ~ (f> for 
^> \XS\. The equations of motion are given by 

4> + 3H<P + m 2 J(f> - (f> ) + X 2 S 2 (f> = 0, (5) 
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S + 3HS + (m 2 s + A 2 2 )S = 0, (6) 

where H is the Hubble parameter. In order to give the large effective mass to enough 
to overcome the Hubble parameter, we assume that S is initially far displaced from the 
origin. Then the effective mass of is given by ~ AS* for |AS| 3> m^. Hence, the 
required condition for to follow the time- dependent minimum adiabatically is |AS| 3> H. 
Under these assumptions, we follow the field dynamics by solving Eqs. (5) and (6) and 
calculate the resultant abundance. 

2.2 The scalar field dynamics 

We investigate the dynamics of the scalar fields after the field S starts to oscillate. For 
simplicity, we assume that the Universe is inflaton oscillation-dominated in the era we 
are interested in and the reheating takes place well after that. 1 Because is initially 
settled on the origin, the mass of S is given by ms, and hence the S begins to oscillate 
at H — ms- The solution of the equation of motion (6) after S starts to oscillate is 

S = A s cos( m s t) with A s = Si — , (7) 

m s 

for H < ms where Si represent the initial amplitude of S, which we regard as a free 
parameter. After S starts to oscillate, its amplitude decreases in proportional to H. 
Therefore, the effective mass of the 0, m^ ff ~ XAs, always exceed the Hubble parameter 
if the following condition is satisfied : A Si > ms- Then, based on the similar arguments 
performed in Refs. [5, 6], the abundance is significantly suppressed by the adiabatic 
suppression mechanism, if the S decay rate (Tg) is so small that S decays after the 
relaxes to the minimum. This condition is rewritten as Ts < = msm^/(XSi). If S 
has a large mass and decays well before BBN, it does not cause significant cosmological 
effects. Then cosmological problems associated with coherent oscillation is avoided. 

Note, however, that there is a subtlety in the model (3). After S starts to oscillate, until 
its amplitude reduces to A s ~ m^/X, the dynamics of is dominated by the interaction 
term A 2 S 2 in the equation of motion (5). Since the universe is assumed to be matter 
dominated, (5) is rewritten as 

X 2 A 2 

u+ — 2~^[1 + cos(2m s t)]u = 0, (8) 

where u = (t/tj)(j) with the initial time £, chosen freely. This equation is mathematically 
equivalent to the Mathieu equation and there exists an exponentially increasing solution. 2 
The behavior of the solution obeys the broad resonance regime discussed in Ref. [10]. The 
number density of u increases exponentially while the adiabatic condition is violated. It 

1 Both <p and S are not assumed to be the inflaton. 

2 We consider the resonant amplification of the homogeneous mode. The finite wave number mode will 
also be amplified, but it does not affect the following estimation of the cf> energy density beyond an 0(1) 
factor. 
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takes place when S crosses the origin, 3 or mst ~ (n + l/2)7r for an integer n, and the 
number density increases by a factor exp(27r/x), where /i is a randomly- changing instability 
parameter and it takes 0.17 on average. 

The stochastic resonance continues as long as the frequency of u is larger than that of 5", 
and ceases at A s ~ m 5 /A = S en d- Note that, as is already mentioned, we need Si > S en d 
in order for the adiabatic suppression to work. If the number of S oscillations during 
Si > As > S en d is huge enough, the resonance continues until the produced energy 
density becomes comparable to that of 5*. This happens at ~ ms/X ~ Send, hence the 
amplitude may be raised up to that of S in the limit of efficient parametric resonance. This 
is shown in Fig. 1(a), where we have solved the equations of motion (5) and (6) and shown 
the time evolutions of and S with following parameters : mj = 1, = 0.01, 0o = 0.1, 
A = 1, Tg = 10~ 4 andSi = 400 in arbitrary units. 

On the other hand, the resonance ceases soon after S starts to oscillate if the duration 
of the resonant oscillation is short. This is the case if values of Si and S en d are rather close 
to each other. In this case, the resulting amplitude of is highly suppressed compared with 
As- Fig. 1(b) shows the time evolutions of and S with m s = 1, = 0.01, O — 0.1, 
A = 1, Ts — 10~ 4 and Si = 50 in arbitrary units. It is clearly seen that the adiabatically 
follows the minimum of the potential without significant oscillations around the minimum 
although some amount of oscillations are induced via parametric resonance. 

The amplitude of is denoted by 0* = £S en d = £,m s /\ at the time A s = S en d = ms/X, 
where £ represents a suppression factor which is determined by the ratio Si/S en d- We 
have found that the suppression is most efficient (£ <C 1) for Si/S ea d ~ 10 and £ = 1 for 
Si/S end > 100. 

Let us briefly summarize necessary conditions for the adiabatic suppression to work. 

• Si > S cn d = rris/\ : Otherwise, the effective mass for is smaller than the Hubble 
parameter and the adiabatic suppression does not work. 

• Si 100S cn d : Otherwise, the resonant amplification of (j) is efficient. Combined 
with the above condition, the suppression is most efficient for Si ~ lOS'end- 

• T s < H ad i = m^Send/ $i '■ Otherwise, S decays and large effective mass for (f> 
disappears before <fi relaxes to the true minimum. 

2.3 The abundance of the oscillating scalar field 

We have seen that the effect of parametric resonance induces the amplitude of = 0* at 
t = t* where £*(~ \Si/m s ) is the cosmic time when S becomes equal to S en d- Now let us 
see how the amplitude of decreases after that. First we consider the dynamics of the 

3 This is the case for a real scalar S. If S is a complex scalar and S 2 in the potential (3) is replaced 
with IS*! 2 , along with a certain magnitude of U(l) violating term oc (S 2 + S* 2 ), the S may rotate in 
the complex plane and not pass through the origin. In this case, there are no such parametric resonant 
processes. 
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(a) Si = 400 



(b) Si = 50 



Figure 1: Time evolutions of XAs (red dotted line), (blue solid line), and H (green 
dashed line) are shown. The small dotted magenta line and the dashed-and-dotted cyan 
line correspond to mj and respectively. Mass-dimensional values are normalized by 
some arbitrary mass scale. We have taken ms = 1, = 0.01, 0o = 0.1, A = 1, and 
= 10~ 4 (the decay rate of S) in both Figures and Si = 400 in Fig. 1(a) and Si = 50 in 
Fig. 1(b). The initial amplitude of is set to be the minimum of the potential (4). 



fields in the era of S en d 3> As 3> m^/X (or ms 3> XAs ^> m$). Since the frequency of S 
(~ ms) is much larger than that of (~ XAs), we can replace the S 2 with its one cycle 
average (S 2 ) = A 2 s /2 which is proportional to t~ 2 . Therefore the equation of motion of 
becomes a simple homogeneous equation and its solution is given by 



t 



-1/2 



COS 



V2XA S . / t 



(9) 



Note that the amplitude of the 0, A^, is damped as A^ oc t~ x l 2 in this era. Next, when 
XA S m^, the bare mass for dominates and the solution of the equation of motion is 
expressed as 

(10) 



i + — cos(m^), 

v 



where B is a constant of integration. This shows that the amplitude of the oscillation of 
around the true minimum decreases as A^ oc t~ l . 

These results are supported by the numerical calculation shown in Fig. 2(a), where 
the time evolutions of and S are shown for 5", = 100 with other parameters taken to be 
the same as those in Fig. 1. From the figure, we can see that the center of the oscillation 
traces the time-dependent minimum of the potential adiabatically. This behavior cannot 
be realized if XSi < ms, in which the minimum changes while is frozen and starts to 
oscillate at H = m§ with initial amplitude 0o, as shown in Fig. 2(b) where we have taken 
Si = 0.1. 

Now we can calculate the abundance of the oscillation. As we have explained we 
assume that S as well as the inflaton decays much later. In this case, the dynamics is 
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(a) Si = 100 (b) S l = 0.1 



Figure 2: Time evolutions of \A$ (dotted red line), H (dashed green line), and (p — (4>) 
(solid blue line) are shown. The dashed-and-dotted cyan line represents m^. Parameters 
are the same as those in Fig. 1 but the initial amplitude of S are Si = 100 in Fig. 2(a) 
and Si = 0.1 in Fig. 2(b). 



qualitably same as that in Fig. 2(a). In order to estimate the <fi abundance, notice that 
the comoving <p number density after the resonance ends is adiabatically invariant. Thus 
we have n^{t) = n^(t*)[a(t#) / a(t)] 3 where n^(t#) = XS en d4>l/2. By noting that the final 
energy density of 4> around the true minimum is given by p^it) = m^n^(t), we obtain the 
following p^ to entropy ratio, 



m s J \MpJ 8 yM, 

where 7 denotes the suppression factor compared with the case in which the adiabatic 
suppression does not occur, and it is given by 

7 ^£W|A 2 (12) 



m s V0| 

Hence the adiabatic suppression is efficient if 7 1. 

We calculated numerically the time evolution of the energy-to-entropy ratio until it is 
fixed at the reheating and the result is shown in Fig. 3(a). We also show the Si dependence 
of the suppression factor 7 in Fig. 3(b). We have taken ms = 1, = 0.01, (fio = 0.1, 
A = 1, and Ts = 10~ 4 . We found that if Si > 200, the resonance becomes efficient and 
£ ~ 1. The resonance is not efficient (i.e. ( C 1) for S% ~ 10, where the suppression 
mechanism works most efficiently. In contrast, for Si < 1, the suppression mechanism no 
longer works. Note that the efficiency of the suppression depends on Si/ S en <i. Therefore, 
for different choice of parameters, different values of Si are favored from the viewpoint of 
the adiabatic suppression mechanism. 
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Figure 3: Time evolutions of p^/s (left) and Si dependence of the suppression factor 7 
(right) are shown. Mass dimensional values are normalized by M = O.OlMp. We have 
taken 171$ = 1, = 0.01, 0o = 0.1, A = 1, and r$ = 10 -4 . In the left figure, we have 
taken Si — 1 (solid red line), Si — 5 (dashed green line), Si = 10 (dotted blue line), 
Si = 100 (small-dotted magenta line), and Si = 200 (dashed-and-dotted cyan line). In 
the right figure, the dashed green line represents the analytical formula given by (12) with 

c 1. 



3 Application to a supersymmetric axion model 

The supersymmetric (SUSY) axion model is introduced to solve the problems of the 
standard model (SM) of particle physics. One of the serious problem is known as the strong 
CP problem. Quantum chromodynamics (QCD) allows the existence of the CP violating 
term in the Lagrangian, but on the other hand, the measurement of the neutron electric 
dipole moment shows that CP must be preserved with high accuracy [11]. This implies 
that the CP violating term must be highly suppressed, and the SM cannot explain this fact 
in a natural way. The most popular solution was proposed by Peccei and Quinn [12]. They 
introduced an additional global U(l) symmetry, called Peccei-Quinn (PQ) symmetry, 
written as {7(1)pq. When the PQ symmetry is broken spontaneously, the axion arises 
as a pseudo Nambu-Goldstone boson [13]. The axion acquires its mass through the 
QCD instanton effect and settles down to the CP preserving minimum of the potential. 
Another problem of the SM is the gauge hierarchy problem. It is naturally solved in the 
framework of supersymmetry (SUSY) [14]. Therefore, the SUSY axion model solves both 
of the problems in the SM. 

In the SUSY axion model, the scalar partner of the axion, saxion, and the fermionic 
superpartner of the axion, axino, are included in the PQ supermultiplet and they take 
significant roles in cosmology [15, 16, 17, 18, 19, 20, 21, 22]. In this section, we focus on 
the dynamics of the saxion and verify that the mechanism shown in the previous section is 
applicable. This was partly mentioned in Ref. [9] in the context of hybrid inflation model 
in the SUSY axion model. We further investigate this issue in a more general form. 
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3.1 The potential of the SUSY axion model 

Let us consider the following superpotential, 

W = K S(VV-f a ), (13) 

where S is a gauge singlet superfield and has a zero PQ charge and *Sf and \l/ are the 
gauge singlet PQ superfields whose PQ charges are +1 and —1, respectively. The PQ 
superfields contain the axion (a), saxion (a) and axino (a). Here f a is the PQ symmetry 
breaking scale and k is dimensionless coupling constant which is assumed to be real and 
positive. The scalar potential derived from the superpotential (13) is 

V = k 2 |M - f a \ 2 + k 2 \S\ 2 (\^\ 2 + |^| 2 ), (14) 

where we used the same symbol for the superfield and its scalar component. The global 
minimum lies at S — and = f 2 . The flat direction along = f 2 is stabilized at 
|^| ~ I 1 !'! ~ f a by including the SUSY breaking effect and correspondingly the saxion, the 
flat direction in the scalar potential of \l/ and obtains a mass of order of the gravitino, 
1713/2. Note that for \S\ > f a , the minimum of the PQ scalar \l/ (\1/) is \l/ = (\I/ = 0), so 
the PQ symmetry is restored. For \S\ < f a on the other hand, the PQ scalars settle on 
the flat direction tyty = f 2 and the PQ symmetry is spontaneously broken. In order to 
apply the adiabatic suppression mechanism to the SUSY axion model, we assume that S 
is initially displaced far from the origin. In order to realize such a condition, for example, 
we add a negative Hubble-induced mass term of S to the potential, which comes from the 
supergravity effect [23], and introduce an additional nonrenormalizable superpotential 

V" en— 1 

W/ NR= 1? - r , (15) 

where X is an additional superfield and M is a cutoff scale and n is a positive integer 
(n > 4). It is achieved by assigning appropriate discrete -R-charges on S and X. Then S 
obtains the potential 

IC|2(n-l) 

V s = -c H H*\S\' + L±^, (16) 

where we assume that X sits at the origin due to the positive large Hubble-induced 
mass. Hence the S tracks the temporal minimum given by \S\ ~ (HM n ~ 3 y/( n - 2 \ which 
is displaced far from the origin. According to the potential (14), the PQ fields ^ and 
^ have the same mass m$ ~ k\S\ soon after the PQ symmetry breaking. Hence the 
expectation values during this period are given by = ( 1 l r ) = f a . The condition that 
the effective mass of the PQ field must be larger than the Hubble parameter is nf a > Hpq 
where Hpq denotes the Hubble parameter at the beginning of S oscillation induced by 
the PQ symmetry breaking. 4 On the other hand, the condition that the energy density 
of the Universe is dominated by the inflaton at the PQ symmetry breaking is written 
as n 2 fa < 3MpHpQ. These conditions are rewritten as f a /M < k < y/3M P /M by 



For n = 4, it is estimated as Hpq ~ for f a <§; kM and Hpq ~ nj a for /„ ^> kM. 
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considering the superpotential (15) with n — 4. Therefore, for a natural value of k, the 
adiabatic suppression is realized. 5 

The dynamics is somewhat complicated. Let us briefly summarize the dynamics before 
going into details described in the following subsections. 

• First, S has a large value due to the potential (16). PQ scalars are stabilized at the 
origin \I> = \I/ = due to the large mass induced by S. S gradually decreases due 
to the Hubble expansion. 

• At H = Hpq, where S ~ f a , the tachyonic instability develops for PQ scalars and 
they fall down to the flat direction tyty = f%. 

• Almost simultaneously, the 5* field begins to feel large mass from the coupling to 

and and then S begins to oscillate around the minimum S = with initial 
amplitude of S ~ f a . 

• Then, the system resembles the toy model in the previous section. The saxion, 
corresponding to the flat direction in the \1/ and \I/ space, obtains a large mass from 
the oscillating scalar S, and it follows the temporal minimum of the potential if the 
following condition is satisfied : nf a ^> Hj>q. 

• The parametric resonance occurs soon after the 5* starts oscillates and it ends within 
one Hubble time. During this short period, however, the resonance effect may 
efficiently amplify the saxion abundance since the oscillation frequency is much 
larger than the Hubble scale. 

• After the resonance ends, the saxion adiabatically follows the minimum until it 
relaxes to the true minimum when 5* reduces to S ~ S^di = m 3 / 2 //t or if ~ H a di = 
HpQ m 3/2/ f a ), if S decays much later. 



3.2 The dynamics of the saxion 
3.2.1 Soon after the PQ symmetry breaking 

First, we consider the PQ field dynamics in some short epoch after the PQ symmetry 
breaking, in which the low energy SUSY breaking terms are assumed to be completely 
negligible : kS ^> m 3 / 2 . We denote the flat direction and the direction perpendicular to 
it as and \l/ 2 , respectively. These directions are given by 

*' = ^5-' *' = (17) 



5 If S dominates the Universe, inflation takes place and ends at S = f a at which tachyonic instability 
along the direction of PQ scalars develops. This is the hybrid inflation model described in Ref. [9]. 
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and the potential of the PQ fields (14) is rewritten as 

v = jf a + K 2 (isi 2 + an 2 + * 2 (\s\ 2 - f 2 a )m\ 2 

K 2 K 2 K 2 (18) 

+ y l*i| 4 + y l* 2 | 4 - y |*i| 2 |* 2 | 2 [l - 2sin 2 (^ + 6 2 )] + Vg, 

where Q x = arg(^i) and 9 2 = arg(# 2 )- Assuming that = ^ 2 = * 2 , and S = S* 
and defining = a/2^i, ip 2 = V%^2, and if = y/2S, the potential is rewritten as 

V = + ^fM - tt) + + + 1^ + - itfty^ + V S . (19) 

From this potential, the equations of motion of respective directions are given by 

(p + 3H<p + y (^ 2 + V 2 > - c H R\ + ^ = 0, (20) 

fa + 3^i + ^ 2 [fl + y - y)^i + y^? = °' (21) 

A + 3i^2 + ^ 2 ( - /a + y - y ) ^2 + y^ 2 3 = 0, (22) 

where we have taken n — 4. Since after the PQ symmetry breaking the PQ scalars 
have (^) = (\E r ) ~ / a , the expectation values of respective directions are (ipi) ~ and 
(ip 2 ) ~ 2/ a . By looking at the mass terms for ipi and ip 2 , it is clear that the ip 2 starts to 
oscillate with initial amplitude of ~ f a around its minimum after the PQ breaking, while 
the flat direction ipi remains at the origin. The approximate solutions of equations (20) 
and (22) are written as 



fa { tj \ . m ~ „ ~ » / 1 



tp-^ljj cos[2Kf a (t - U)], i) 2 ~ 2f a + 2f a lj\ cos[2Kf a (t - *,-)]■ (23) 

Looking at the equation of motion of the flat direction (21), there exists a interaction 
term which may induce the parametric resonance, as shown in the previous section. Since 
ip ~ f a soon after the PQ breaking, the frequency of the oscillation %p\ is almost equal 
to that of tp2i so the parametric resonance takes place within a short epoch after the 
PQ breaking. However, the frequency itself is much larger than the Hubble scale at that 
epoch (nfa ^> H?q), so there may be many oscillations during this short period. Thus 
we cannot necessarily neglect the effect of parametric resonance. 

To verify this, we have numerically solved the set of equations of motion (20)-(22). 
The numerical results are shown in Fig. 4, where time evolutions of Ktp (small dotted 
magenta line), ip 2 (dashed red line), and ipi (solid blue line) are shown. We have taken 
ip = 2, H = 6 x 10~ 3 , and ip\ = ip 2 = 10~ 20 as initial values and k — 1 for Fig. 4(a) 
and k = 0.01 for Fig. 4(b). In fact, the initial value of ip\ is not reflected unless it is 
too large because the initial amplitude of ip\ is determined by the tiny deviation from 
^ = ^ = f a at the PQ breaking, which is induced by the gravitino mass terms. The 
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(a) 



(b) k = 0.01 



Figure 4: Time evolutions of K(p (small dotted magenta line), ?/> 2 (dashed red line), and 
ipi (solid blue line) are shown. Mass dimension is normalized by the PQ scale f a . We 
have taken (p = 2, H = 6 x 10~ 3 , and ip\ = ip2 — 10~ 20 as initial values and 1713/2 = 10 -8 , 
n — 4, ch — 1, M — 5 x 10 4 for both Figures, and k = 1 for Fig. 4(a) and k = 0.01 for 
Fig. 4(b). 



resonance is efficient only if the Hubble parameter is much less than nf a at the beginning 
of the oscillation of <p and ip2- We investigated it numerically and concluded that if 
Hpq < 10~ 2 fi;/ a , the amplitude of the flat direction is raised up to f a . In the case of 
Hj>q ~ 10 _1 /t/ a , the resonant amplification of ip\ becomes inefficient much sooner due to 
the expansion of the Universe, and the amplitude of the flat direction is highly suppressed 
as shown in Fig. 4(b). The saxion is identified as a oscillation along the flat direction, so 
the initial amplitude of the saxion can be expressed as cr, = £/ a with £ < 1. In particular, 
in the case of Fig. 4(b) (Kf a /HpQ ~ 10 at the beginning of S oscillation and k = 0.01), 
we found £ ~ 10 -8 . 



3.2.2 Well after the PQ symmetry breaking 

As the Universe is cooled, the low scale SUSY breaking terms become dominant. It 
happens at S < S^di = ^13/2/ K when the Hubble parameter becomes smaller than H^i = 
HpQfn 3 /2/(Kfa)- Then, the dominant contribution to the potential of the PQ fields along 
the flat direction is 

\/ = fi: 2 |5| 2 (|^| 2 + |^| 2 ) + m 2 /2 (c 1 |^| 2 + c 2 |^| 2 ) with M = / a 2 , (24) 

where c\ and c 2 are assumed to be real and positive constants of order unities and we 
neglect the contribution of (16). The time-dependent minimum for \l/ is derived as 

//€ 2 |S| 2 + c 2 m 2 \V4 

(*)=( 2 LL f fa- (25) 
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The equations of motion are given by 

S + 3HS + k 2 {\V\ 2 + \V\ 2 )S = 0, (26) 
* + 3H^ + [k 2 (\S\ 2 + |^| 2 ) + c im l /2 ]y - K 2 f 2 ^* = 0, (27) 
§ + 377^ + [k 2 (|S| 2 + |^| 2 ) + c 2 m 2 /2 ]^ - k 2 / 2 ^* = 0. (28) 

Let us redefine the PQ field along the flat direction as a = 2[Re(\I / ) — (\&)] and call it 
saxion. Using this definition, well after the PQ fields approach the minimum, the potential 
of the saxion is written as V = A(k 2 \S\ 2 + cim 2 ^ 2 )a 2 , so the equation of motion of the 
saxion becomes 

a + 3H& + 8(k 2 \S\ 2 + c x m\ l2 )a + 2(tf ) + 6H(V) = 0. (29) 

The behavior of the solution of this equation of motion is similar to that of the toy 
model discussed in the previous section. While kS ^> m 3 / 2 , the amplitude of the saxion 
decays like t _1//2 as long as the Hubble parameter does not exceed the effective mass of 
the PQ fields. After the saxion mass (7713/2) becomes significant (777.3/2 ^> kS), the saxion 
amplitude decreases as t -1 . The numerical calculations are shown in Fig. 5, where time 
evolutions of scalar fields are plotted. Similarly to the toy model, the saxion follows its 
minimum adiabatically even after the temporal minimum shifts to the true minimum 
at kS ~ m.3/2. Thus the resulting amplitude of the saxion is significantly suppressed 
compared with the ordinary case. 

Let us briefly summarize necessary conditions for the adiabatic suppression to work. 

• K fa ^ Hpq : Otherwise, the effective mass for saxion is smaller than the Hubble 
parameter and the adiabatic suppression does not work. 

• K fa ^ lOOilpQ : Otherwise, the resonant amplification of saxion is efficient. Com- 
bined with the above condition, the suppression is most efficient for nf a ~ lOHpq. 

• < if a di = H P Qm 3 / 2 /(Kf a ) : Otherwise, S decays and large effective mass for the 
saxion disappears before the saxion relaxes to the true minimum. 

In the next subsection, we calculate the abundance of the saxion and see how much 
the saxion coherent oscillation survives. 

3.3 The abundance of the saxion 

As shown in the previous work [21], the abundance of the saxion is strongly constrained 
from cosmological observations. If we succeed in reducing the abundance of the saxion, 
the cosmologically viable model parameter space can be extended, and it is very important 
when one tries to build a SUSY axion model. Now, using the setup we developed in the 
previous subsections, we calculate the abundance of the saxion and see that the adiabatic 
suppression significantly reduces the final abundance. 
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Figure 5: The time evolutions of \I/ — (displacement from the temporal minimum, 
thin solid cyan line), \I/ — ^ f (displacement from the eventual minimum, thick solid red 
line), and nAs (dotted blue line). The Hubble parameter and the gravitino mass are 
represented as dashed green line and dashed-and-dotted magenta line respectively. Mass 
scale is normalized by f a . We have taken S — 1, ^ = 0.1 + {^), and H = 0.006 as initial 
values and k = 1 and m 3 / 2 = 0.01 




The abundance of the saxion depends on the decay rate of S and decay rate of the 
inflaton. Here we assume that the decay of S and inflaton occurs after the saxion relaxes 
to the true minimum, i.e., T s , Tj < if adi = H P Qm 3 / 2 / (nf a )- In particular, the former 
condition is necessary for the adiabatic suppression works successfully. These conditions 
are written as 

m 3/2 > 7^r 5 and m 3/2 > Tp^Tj. (30) 

-HPQ -HpQ 

The latter condition is rewritten as 

( Kfa \ ( T R 



m 3/2 > 10- b GeV Ja . (31) 

d/ ~ \100H PQ J Vl0 5 GeVy V ; 

Under these assumptions, the energy-to-entropy ratio of the saxion is fixed at the 
reheating (H ~ Tj). By noting that the comoving saxion number density after the 
resonance ends is adiabatically invariant, the number density is calculated from n a (H) = 
n(T (H PQ )[a(H PQ )/a(H)] 3 for H < H PQ , where n a (H PQ ) = nf a af/2 = Kf a (£f a ) 2 /2. The 
energy-to-entropy ratio is then estimated as 
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Figure 6: Time evolutions of p a js and k dependence of the suppression factor 7 are 
shown. Mass dimensional values are normalized by /„. We have taken /„ = 10 12 GeV, 
m 3/2 = 10" 5 , r s = r> = 10" 6 , A = 10~ 5 , n = 4,c H = l, and M = 10 3 in both figures. In 



the left figure, we have taken k — 1 (solid red line), k = 0.1 (dashed green line), k 
(dotted blue line), k = 0.001 (small-dotted magenta line). 



0.01 



where 7 is a suppression factor defined by 

7 = e (^) 2 ^Il = m-Y'(^hX(°±) ( Wl \ GeY ) |. (33) 



HpqJ Ufa \ W0H PQ J \ k I V fa J \l GeV 

Here A represents the difference between the high-energy and low-energy minima, which 
is written as A ~ 1 1 — (C2/ ci) 1//4 |, and this may be order unity in general. If there was no 
adiabatic suppression and the parametric resonance, we would have 7 = 1. Therefore sig- 
nificant suppression is certainly realized unless nf a is too large compared with the Hubble 
parameter at the PQ breaking. Note that the most effective suppression takes place if 
Kf a /HpQ ~ 10 because the parametric resonance ends soon and the saxion amplitude is 
not amplified much. 

We calculated numerically the evolutions of the energy-to-entropy ratio and K-dependence 
of the suppression factor 7, and the results are shown in Figure 6. In this calculation, 
mass dimensional values are normalized by f a and in order make numerical calculation 
easy we have taken /„ = 10 12 GeV, m 3/2 = 10" 5 , T 5 = T 7 = 10~ 6 , A = 10" 5 , n = 4, 
c H = 1, and M = 10 3 in Figure 6. We found that nf a /H PQ ~ 10 3 , 10 2 , 10, 1 for 
k — 1, 0.1, 0.01, 0.001 respectively in this calculation. We can see that the most efficient 
suppression takes place for k ~ 0.01 and the resonant amplification becomes efficient for 
n > 0.1. For k < 0.001, the adiabatic suppression no longer works and 7 approaches 
unity. Note that 7 becomes much smaller than unity in realistic models even if the res- 
onant amplification takes place most efficiently, because the realistic gravitino mass is 
many orders of magnitude smaller than that we use in our numerical calculation. 

The saxion is also produced from thermal bath and the contribution to the saxion 
abundance should not be neglected. The temperature at which the saxion decouples from 
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the thermal bath is estimated as [15] 



10 12 GeV 

The abundance of the thermal saxion is then estimated as 

(TP) 



~10- 3 Gevf^J for T R >T D , (35) 



and 

These predictions are summarized in Fig. 7. In our model, the contribution of the saxion 
coherent oscillation to the total abundance can be much smaller than that from the ther- 
mal scattering. In particular, considering the case of the small £, in which Kf a /Hpq ~ 10 
and k = 0.01, the abundance of the coherent oscillation becomes much smaller and com- 
pletely negligible compared to the thermal saxion. Looking at Fig. 7(b), it is remark- 
able that the cosmological constraints for the saxion is much weakened than previously 
thought [21]. 

Note that since the PQ symmetry is broken after inflation, cosmic strings are formed 
at the PQ phase transition. Moreover, at the QCD phase transition, domain walls ap- 
pear [24]. In order to avoid the cosmological domain wall problem, the color anomaly 
number must be equal to one. In this case domain walls are bounded by strings, and 
walls as well as strings disappear due to the wall tension [25, 26, 27, 28, 29]. The PQ 
scale, f a , cannot be larger than ~ 10 12 GeV in this case since axions produced through 
coherent oscillation and emission from the strings and walls contribute to the dark matter 
density [30, 31]. 



4 Conclusions 

We revisited the problems induced by the scalar coherent oscillation such as the moduli 
or saxion. In the ordinary consideration, the scalar field initially displaced far from the 
origin begins to oscillate when the Hubble parameter becomes equal to the mass of the 
scalar field, which leads to the large abundance and often becomes harmful for cosmology. 

In this paper, we have proposed a model in which the amplitude of oscillating scalar 
fields is significantly suppressed by using the idea of adiabatic suppression. This mecha- 
nism works if the effective mass of the scalar field is much lager than the Hubble parameter. 
We found that such a situation is realized by coupling the scalar field with another scalar 
field, whose initial amplitude is also far displaced from the origin. We have investigated 
these scaler field dynamics and seen that the oscillating scalar field follows the temporal 
minimum of the potential adiabatically without inducing a large oscillation amplitude 
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(a) /„ = 10 12 GeV (b) /„ = 10 14 GeV 

Figure 7: The mass dependence of the energy-to-entropy ratio of the saxion predicted 
by our analysis. The solid blue lines represent the our results for the saxion coherent 
oscillation, the dashed green lines represent ordinary analysis of the saxion coherent os- 
cillation without the adiabatic suppression (7 = 1), and the red dotted lines represent 
the thermally-produced saxion. Thick and thin lines correspond to the reheating temper- 
ature Tr = 10 5 GeV and T R — 1 GeV respectively. In this calculation, we have taken 
KfjHpQ = 100, k = 0.1, f = 1, and f a = 10 12 (10 14 ) GeV for left (right) figure. The 
dashed-and-dotted blue line represents the most suppressed case of coherent oscillation, 
in which we have taken nf a /HpQ = 10, k — 0.01, £ = 10~ 8 , and Tr = 10 5 GeV. 

depending on model parameters. In the model we have analyzed, we have encountered 
a parametric resonance for amplifying the oscillation amplitude, but still the abundance 
can be significantly reduced. Although the resulting abundance depends on the model 
parameters, a broad parameter region is available for successful adiabatic suppression. 

We have applied the adiabatic suppression mechanism to the SUSY axion model. We 
followed the dynamics of the PQ scalar fields and calculated the resulting abundance of 
the saxion. According to the result, the cosmo logical constraint for the saxion is greatly 
relaxed particularly for large PQ breaking scale without relying on any additional entropy 
production. Although we have focused on the SUSY axion model as a concrete example, 
we believe that the idea has a broad applicability. A similar mechanism may be used 
for solving the cosmological problems associated with scalar field oscillation appearing in 
many models beyond SM. 
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